Run the scripts in the following order on Comsol livelink with Matlab:

create_model.m: create the comsol model (geometry, physics, mesh) (simulation time is very long; recommended to use a cluster)
solve_model.m: solve the comsol model (it will require a large amount of memory (greater than 500 GB), hence one needs to use a high memory cluster to run this script)
plot_mode.m: plot the mode profile of the output modes (need to first run the above two scripts)